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Introduction 


The  long-term  goal  of  this  research  is  to  develop  a  non-invasive  technique  for  clinical  management  of  TSC- 
associated  renal  angiomyolipomas.  TSC  is  a  genetic  tumor  predisposition  syndrome  characterized  by  the 
growth  of  lesions  in  multiple  organ  systems.  Approximately  80%  of  TSC  patients  develop  renal 
angiomyolipomas,  a  type  of  lesion  composed  of  variable  amounts  of  fat,  smooth  muscle,  and  vascular  tissue. 
Renal  angiomyolipomas  are  often  benign  and  present  with  multiple  lesions  in  each  kidney.  Patients  with  renal 
angiomyolipomas  may  experience  discomfort,  flank  pain,  hydronephrosis,  hematuria,  and  hypertension.  These 
lesions  can  also  lead  to  acute  hemorrhaging  or  chronic  loss  of  renal  function.  In  this  research,  we  will  use  High 
Intensity  Focused  Ultrasound  (HIFU)  to  ablate  tumors  and  magnetic  resonance  (MR)  imaging  to  monitor 
ablation.  MR-guided  HIFU  enables  “surgical  procedures”  to  be  performed  deep  within  the  body  without 
incisions  or  punctures,  providing  a  risk-free  therapeutic  approach  to  managing  benign  or  malignant  lesions. 

The  proposed  project  aims  to  investigate  whether  MR-guided  HIFU  offers  the  ability  to  ablate  renal 
angiomyolipomas  in  a  mouse  tumor  model.  We  expect  that  the  research  experience  gained  from  this  pre- 
clinical  study  will  benefit  our  future  study  on  non-invasive  thermal  therapy  of  renal  angiomyolipomas  in  human 
with  MR-guided  HIFU.  In  the  proposed  work,  we  had  planned  to  focus  on  the  development  of  an  MR-guided 
HIFU  experimental  system  for  thermal  ablation  in  mouse  kidneys  in  the  first  year  and  perform  animal 
experiments  using  the  technique  developed  in  the  second  year. 

The  physical  mechanisms  underlying  HIFU  is  that  a  HIFU  transducer  constructed  with  a  concave  shape  and/or 
multiple  elements  has  the  ability  to  focus  acoustic  energy  into  a  target  volume  having  a  diameter  of  a  few 
millimeters.  The  focused  acoustic  energy  induces  a  rapid  rise  in  temperature  (e.g.  70°C  to  100°C),  resulting  in 
thermal  necrosis  of  tissues  in  the  target  volume.  Although  HIFU  offers  the  capability  of  thermal  ablation,  non- 
invasive  thermal  therapy  is  possible  only  if  the  focal  spot  of  HIFU  can  be  controlled  within  the  body  using  the 
feedback  information  provided  by  medical  imaging  guidance.  MR  is  superior  to  other  imaging  modalities 
because  it  provides  both  excellent  soft-tissue  visualization  and  the  ability  to  monitor  thermal  delivery 
(temperature  mapping).  However,  an  intrinsic  drawback  of  MR  is  the  slow  data  acquisition  speed  limited  by 
long  relaxation  time.  As  a  result,  we  had  expected  to  meet  challenges  in  imaging  mouse  kidneys  with 
respiratory  movement. 

Our  first  year's  work  led  to  the  development  of  a  small  animal  HIFU  system  with  feedback  control.  In  several  in- 
vitro  studies  without  motion  concerns,  we  have  demonstrated  that  the  developed  HIFU  system  offers  the  ability 
to  focus  acoustic  energy  within  a  small  spot  in  a  diameter  of  a  few  millimeters.  During  the  second  year,  we 
proceeded  with  in-vivo  studies  and  found  that  motion  may  be  destructive  to  HIFU  thermal  treatment.  Although 
parallel  imaging  has  been  proposed  and  used  to  address  this  issue  in  this  work,  it  was  found  that  the 
performance  of  parallel  imaging  is  not  satisfactory  because  coil  array  configuration  is  limited  by  small  mouse 
anatomy.  This  challenge  has  delayed  our  animal  experiment  and  we  are  actively  seeking  a  different  solution  to 
resolving  the  problem.  The  proposed  project  has  been  extended  until  2014  and  a  different  research  plan  will  be 
used  to  advance  our  research. 

This  report  is  a  summary  of  our  second  year's  work  on  the  proposed  research.  The  summary  will  focus  on  the 
technical  challenges  we  met  and  provide  details  about  why  this  challenge  cannot  be  overcome  using  the 
proposed  techniques.  In  addition,  we  will  propose  our  alternative  plan  to  conduct  the  proposed  animal 
experiments  in  the  third  year.  In  this  plan,  two  methods  will  be  used  to  overcome  the  technical  challenge  we 
found  in  the  second  year's  work. 
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Body 


In  the  second  year,  our  effort  was  focused  on  the  proposed  Task  3.  We  found  that  the  proposed  parallel 
imaging  approach  could  not  meet  the  requirements  for  feedback  control  in  in-vivo  studies  due  to  the  limitation 
posed  by  small-animal  anatomy  and  acoustic  power  transmission. 

Proposed  Task  3:  Optimization  of  the  MR-guided  HIFU  system  for  thermal  ablation  of  renal 
angiomyolipomas  in  a  mouse  orthotopic  xenograft  model  (Stage  2:  Months  13-18). 

The  real-time  MR  guidance  of  renal  targets  will  be  optimized  in  a  free-breathing  mouse  orthotopic  xenograft 
model  developed  in  Task  2.  The  adaptive  control  system  will  be  optimized  for  maximizing  the  effectiveness  and 
safety  of  HIFU  thermal  ablation.  10  mice  will  be  used.  This  in-vivo  test  is  the  preparation  phase  for  the 
following  systematic  evaluation  of  MR-guided  HIFU  therapy  in  task  4. 


Research  accomplishment:  Figure  1  shows  the  experimental  system  we  developed  for  thermal  ablation  in 
phantoms  and  animals  during  the  first  year  of  this  project.  This  setup  allows  the  HIFU  console  to  synchronize 
MRI  scanning  when  running  HIFU  thermal  ablation.  The  acoustic  power  delivery  from  the  HIFU  transducer  can 
be  dynamically  updated  by  the  HIFU  console  based  on  MRI  information,  providing  a  feedback  control  for  the 
thermal  delivery  deep  within  the  body.  During  second  year,  we  used  this  system  to  investigate  whether  HIFU 
focal  lesions  may  be  formed  deep  inside  the  mouse  body.  In  our  initial  attempts,  healthy  mice  were  used. 


(a) 


MRI  reception  & 
image  formation 


Thermal  Dose  control  & 
therapy  optimization 


Figure  1.  A  feedback  control  system  for  preclinical  thermal  therapy  was  established  using  MR-guided  HIFU 
in  the  CCHMC  IRC.  The  feedback  control  loop  is  formed  using  a  3.0  Tesla  32-channel  Philips  Achieva  MR 
scanner  (Philips  Healthcare,  Best,  the  Netherlands)  and  a  Philips  small-animal  HIFU  system  (Philips 
Healthcare,  Vantaa,  Finland).  This  system  uses  an  eight-channel  3  MHz  sectored  transducer  (IMASONIC, 
Voray  sur  I'Ognon,  France)  for  acoustic  transmission,  (a)  Experimental  setup,  (b)  An  MRI  coil  was  built  for 
mouse  imaging  on  Philips  3T  scanner,  (c)  A  mechanic  stage  was  constructed  for  holding/stabilizing  the 
mouse  and  the  coil  within  the  MRI  scanner.  Inside  the  stage,  water  will  be  used  as  the  interface  between  the 
HIFU  transducer  and  the  animal.  MRI  console,  (d)  Placement  of  transducer  and  mouse  inside  the  MRI 
scanner  using  the  mechanic  stage  in  (c). 


2a.  Optimization  of  real-time  MR  guidance  (Months  13-16). 

Parallel  imaging  will  be  optimized  to  improve  imaging  speed  and  minimize  motion  artifacts  in  free-breathing 
mice.  Specifically,  the  optimization  will  be  performed  to  improve  four  major  imaging  methods,  Ti  weighted 
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imaging,  T2  weighted  imaging,  stiffness  weighted  imaging,  and  phase  imaging.  It  will  be  demonstrated  that 
these  imaging  methods  can  provide  accurate  and  real-time  information  about  HIFU  therapy  delivery. 


Research  accomplishment:  In  this  work,  we  found  that  the  proposed  parallel  imaging  technique  cannot 
provide  good  performance  due  to  coil  array  limitation.  Figure  2  shows  a  HIFU  ablation  result  using  the 
experimental  system  (Figure  1)  and  a  healthy  animal.  It  was  found  that  HIFU  energy  was  successfully 
delivered  into  the  animal  body  while  temperature  information  was  not  able  to  be  collected  in  real  time  using 
parallel  imaging  techniques.  Our  further  investigation  shows  that  the  coil  array  we  developed  (Figure  1b) 
cannot  provide  good  parallel  imaging  performance  because  the  coil  elements  have  to  be  larger  than  the  mouse 
anatomy  and  positioned  in  a  certain  distance  to  the  mouse.  This  configuration  was  practically  limited  by  two 
factors:  First,  acoustic  power  transmission  requires  an  open  space  to  the  abdominal  region  of  the  mouse. 
Since  water  has  to  be  used  to  provide  acoustic  interface  between  the  transducer  and  the  mouse  body,  coil 
elements  cannot  be  positioned  close  to  the  imaging  target  zone.  Second,  the  coil  elements  have  to  large 
enough  to  achieve  certain  RF  penetration  due  to  the  positioning  limitation.  As  a  result,  spatial  variation  of  multi¬ 
channel  coil  sensitivity  was  low,  which  significantly  reduces  the  spatial  encoding  of  coil  sensitivity  within  the 
mouse  body.  Without  sufficient  imaging  speed,  dynamic  imaging  provides  low  image  quality  and  temperature 
information  (phase  change  of  MR  images)  was  distorted  by  imaging  artifacts. 


ir  '  I 
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T1  image  before  ablation 

T 1  image  after  ablation 

Figure  2.  Mouse  HIFU  ablation.  Using  the  experimental  system  in  Figure  1,  HIFU  ablation  was  performed  in 
a  healthy  mouse.  T1 -weighted  MR  images  were  collected  to  determine  whether  a  HIFU  lesion  was  formed 
in  the  target  zone.  By  comparing  T1 -weighted  images  before  and  after  ablation,  the  contrast  difference  was 
found,  implying  HIFU  lesion  was  produced  deep  inside  the  animal  body.  However,  dynamic  MR  images 
were  not  able  to  provide  accurate  temperature  information  for  real-time  feedback  control. 


2b.  Optimization  of  the  adaptive  control  system  for  HIFU  (Months  17-18). 

In  free-breathing  mice,  the  adaptive  control  system  will  be  optimized  to  improve  the  localization  and  power  of 
HIFU  focal  spot  dynamically  based  on  real-time  MR  guidance.  It  will  be  demonstrated  that  HIFU  lesions  can  be 
effectively  created  irrespective  of  respiratory  motion.  The  safety  limit  for  thermal  dose  in  HIFU  ablation  will  also 
be  determined.  This  safety  limit  will  be  used  to  determine  the  thermal  dose  for  MR-guided  HIFU  therapy  in  the 
following  task  4. 


We  have  mathematically  investigated  the  model  for  the  implementation  of  adaptive  control  in  HIFU  treatment 
(See  appendix  1).  In  this  model,  HIFU  was  considered  as  an  inverse  problem.  Our  collaborator  developed  a 
simulated  model  to  test  the  possibility  of  adaptive  control  in  real  time.  We  are  working  on  an  approach  to  use 
the  simulated  model  to  predict  HIFU  treatment  in  mouse. 


Alternative  approaches  to  animal  experiments:  To  overcome  the  challenge  we  found  in  the  second  year's 
work  (Task  3),  we  have  launched  two  efforts  to  seek  a  solution. 

First,  we  developed  research  collaboration  with  Dr.  Qiming  Zhang's  group  in  Pennsylvania  State  University. 
This  group  is  working  on  a  new  type  of  Magnetoelectric  Sensors  for  MR  imaging.  These  sensors  have  the 
potential  to  provide  imaging  performance  similar  to  RF  coil  array.  However,  their  physical  dimensions  are  in 
micro-meter  (Figure  3)  and  they  can  be  positioned  as  close  to  the  mouse  body  as  possible  without  interaction 
on  HIFU  acoustic  transmission.  In  addition,  since  these  sensors  detect  MR  signals  from  sensing  the  mechanic 
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force  generated  by  small  magnetic  field,  no  RF  coupling  is  generated  between  different  sensors.  As  a  result, 
the  channel  number  for  MR  imaging  is  unlimited.  This  is  a  significant  advantage  over  conventional  coil  array 
technology.  We  are  expecting  to  have  the  sensors  ready  for  a  test  in  January  2014.  Once  we  have  a  set  of 
sensors  for  our  HIFU  system,  we  will  be  able  to  continue  our  experiments. 


Our  second  choice  is  to  use  the 
human  HIFU  system  being  installed 
at  our  institute.  We  purchased  a  1 .5 
Tesla  Ingenia  MR  imaging  system 
with  human  HIFU  from  Philips 
Healthcare.  Currently,  the  system 
was  already  delivered  and  is  being 
installed  in  Imaging  Research 
Center  (Figure  4).  This  human  HIFU 
system  has  256  acoustic  elements 
and  offers  the  capability  of  moving 
the  transducers  in  an  arbitrary 
direction.  More  importantly,  the  commercial  HIFU  system  allows  the  use  of  digital  coil  arrays  provided  by 
Ingenia  system.  As  a  result,  we  expect  that  this  new  system  provides  much  better  performance  than  the  small 
animal  HIFU  system  we  are  currently  using.  Once  the  system  is  ready,  we  planned  to  continue  our  animal 
experiments. 

The  final  timeline  of  our  animal  experiments  will  be  dependent  on  the  progress  of  our  two  efforts  described 
above.  It  is  expected  that  further  experiments  on  mouse  will  be  conducted  in  January  2014. 


Key  Research  Accomplishments 

1 .  We  found  a  technical  challenge  that  arises  from  the  coil  array  limitation  posed  by  small  animal  anatomy 
and  acoustic  power  transmission. 

2.  We  launched  two  independent  efforts  to  address  the  technical  challenge  met  in  the  proposed  Task  3. 

a.  Magnetoelectric  sensors  are  being  developed  through  the  research  collaboration  with 
Pennsylvania  State  University. 

b.  A  new  human  HIFU  system  is  being  installed.  We  are  planning  to  transfer  the  animal 
experiments  to  this  system. 

3.  A  mathematical  model  was  developed  to  investigate  adaptive  control  algorithms  in  mouse  HIFU 
ablation. 


Reportable  Outcomes 


1.  We  have  collaboratively  developed  a  mathematical  model  with  Dr.  Donald  French's  group  in  the 
Mathematical  Department  at  the  University  of  Cincinnati  (Appendix  1).  This  model  provides  a  tool  to 
investigate  adaptive  control  problems  in  mouse  HIFU  treatment. 

2.  We  have  developed  our  collaboration  Dr.  Qiming  Zhang's  group.  We  are  collaboratively  developing  a 
new  type  of  magnetoelectric  sensors  that  may  be  used  to  replace  coil  array  in  our  current  HIFU  system 
for  mouse  ablation. 


sensor  we  are  developing  MR  imaging  system  and  human  HIFU  machine  at 
or  MR  imaging.  the  Cincinati  Children's  Hospital  Medical  Center. 


Conclusion 

In  summary,  we  have  found  new  challenges  in  this  project.  From  our  preliminary  investigation  using  in-vivo 
studies,  we  understood  the  physical  mechanisms  underlying  the  technical  challenge  and  launched  two 
independent  efforts  to  address  this  technical  challenge.  We  expect  to  continue  our  animal  experiments  at  the 
beginning  of  2014. 
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1.  Introduction: 

High  intensity  focused  ultrasound  (HIFU)  is  applied  to  a  region  like  the  one  in  figure  1A  which  we  label  O 
and  consists  of  tissue,  skin,  fat  and  blood  in  order  to  ablate  cancerous  tumors.  Magnetic  resonance  imaging 
(MRI),  figure  IB,  is  used  to  measure  tissue  temperatures  and  then  estimate,  with  further  mathematical 
modeling  and  scientific  computing,  various  material  parameters  (e.g.  conductivities,  perfusions  etc.)  in  the 
region  and  monitor  the  tumor  remediation  process.  In  this  draft  we  study  the  mathematical  model  and 
various  estimation  procedures. 

A.  B. 


Figure  1:  Expermental  apparatus  and  diagram  for  fMRI  assisted  HIFU  analysis.  (A)  Example  problem  with 
Ultrasound  transducer  on  the  left  and  potential  treatment  region  on  the  right.  The  diameter  of  the  whole  ’’liver” 
would  be  around  10  cm.  (B)  Magnetic  resonance  imaging  machine  (MRI). 

Following  [HHMK],  we  use  the  Helmholtz  partial  differential  equation  to  model  the  pressure  waves  arising 
from  the  HIFU  transducers; 

/ 1  \  fc2 

V  •  (  -Vp  )  H - p  =  0  in  H.  (1) 

\P  J  P 

Here  k  =  2i xf/c  +  ia,  p  is  a  positive  density  function,  a  is  an  absorption  coefficient,  c  is  sound  speed  and  / 
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is  a  frequency.  We  append  boundary  conditions 


dp 

dv 


i2nfpvn  on  IV  and 


on  dYl  —  IV. 


(2) 


The  transducer  velocity  vn  is  provided  on  boundary  portion  Yu •  The  temperature  T  in  the  liver/tumor 
which  is  induced  by  heating  from  the  pressure  waves  is  read  by  the  MRI  and  evolves  via 


P°T~di  =  V  ’  (/"VT)  “  ^ T  ~  Ta)  +  Q  where  P  =  “bCb-  (3) 

ujb  is  blood  flow  rate;  Ct  &  Cb  are  specific  heat  capacities  of  the  tissue  and  blood,  respectively.  The 
/3(T  —  Ta)  term  models  the  averaged  effects  of  blood  vessels.  Also,  ft  is  a  conductivity  coefficient,  Ta  is  an 
ambient  temperature  and  Q  =  a\p\2/ ( cp ).  Boundary  and  initial  conditions  for  T  =  T(x,  t)  are  added  as  well; 


T(-,  0)  =  Ta  and  T  =  TA  on  dtt. 


(4) 


In  tumor  treatment,  the  values  of  the  absorption,  heat  conductivity,  pressure  density  and  perfusion,  that 
is,  <a,  ft,  p  and  /3  are  usually  unknown.  These  parameters  are  typically  piecewise  constants  dependent  on 
variations  between  air,  skin,  tissue,  blood  and  tumor.  The  main  goal  of  this  work  is  to  estimate  these 
parameters  given  data,  measured  by  the  MRI,  on  temperature  T. 

The  ultrasound  is  crucial  here  since  variations  in  temperature  are  important  in  the  identification  of  the 
unknown  parameters.  These  ultrasound  excitations  cause  temperature  changes  of  3-5  degrees  C.  In  our  lab 
temperature  can  be  measured  on  a  grid  with  0.5  mm  subinterval  widths  (assuming  the  area  being  treated 
is  roughly  10  mm)  and  a  5  second  time  step.  We  believe  it  may  also  be  possible  to  read  pressures  using 
techniques  discussed  in  [GFMME]  or  [LPADSSMBFT].  Of  course,  if  both  temperature  and  pressure  readings 
are  desired,  the  measurements  will  lie  on  coarser  grids. 

The  most  recent  and  related  work,  that  we  are  currently  aware  of,  is  [HHMK].  In  this  paper  only 
temperature  readings  are  available  as  data  and  they  assume  that  a  and  p  are  known.  Their  analysis  area  is 
is  segmented  into  three  zones  so,  overall,  there  are  six  unknowns  constants. 


2.  Focusing  in  Helmholtz  Model: 


The  Helmholtz  partial  differential  equation  is  used  in  [HHMK]  to  model  the  acoustic  pressure  waves  from 
an  HIFU  transducer: 


V- 


/ 1  \  k2 

-Vp  H - p  =  0 

\P  J  P 


in  Cl. 


(5) 


Here  k  =  2i rf/c  +  icq  p  is  a  positive  density  function,  a  >  0  is  an  absorption  coefficient,  c  >  0  is  speed  of 
sound  and  /  is  frequency  of  the  wave  field.  We  append  boundary  conditions 


dp 

dv 


i2irfpv  on  Yu 


and 


on  dYl  —  Yu- 


(6) 


The  velocity  v  is  provided  by  the  ultrasound  transducer  on  boundary  portion  Yu.  In  [HHMK],  the  transducer 
is  pulsed  and,  given  an  array  of  transducers  may  have  many  v  values. 

To  simplify  and  make  some  first  steps  toward  understanding  the  focussing  that  tends  to  occur  in  homo¬ 
geneous  media  we  assume  that  p  is  constant,  let  u)  =  2i r/,  require  a  =  0  and  assume  the  problem  is  posed 
on  a  disk  with  angular  symmetry  and  radius  R,  so  (5)  becomes,  find  p  =  p(r)  for  r  E  [0,  R]  so 


r2p"  +  rp'  +  k2r2p  =  0  with  p'(R)  =  ipvk 
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and  we  assume  \imr^o  (rp'  (r))  =  0  as  a  disc  center  condition.  This  will  have  a  Bessel  function  solution;  to 
derive  this,  we  change  variables.  Let 

r  —  x/k  and  P(x)=p(x/k)  (So  p(r)  =  P(kr)) 


and  find 

x2P"  +  xP'  +  x2P  =  0  and  P'(kR)  =  icpv 
It  is  well  known  (See,  for  instance,  [BC])  that  the  solution  to  the  differential  equation  is 


P(x)  =  AJ0(x)  P  BY0(x). 


But,  since,  Yq  — >  ±00  as  x  — >  0  we  typically  take  B  =  0.  Also,  since  it  is  known  that  Jq  =  —J\  and  thus  we 
can  determine  a  formula  for  A; 


So, 


— AJi(kR )  =  icpv  => 


icpv 

J^kR)' 


=>  «r)—jmMkr) 

From  the  identity  (Again,  see  [BC]) 

2  r/2 

Ji(z )  =  —  /  sin(2:  sin(0))  sin(0)  dcf) 

7T  Jo 


I  believe  we  can  make  a  reasonable  argument  that  J\(kR)  ~  1/k  for  k  large  (I’m  getting  abit  on  a  limb  with 
this  one  . . .)  and  thus  have  \A\  ~  vk.  From  the  picture  below,  we  now  would  have  that,  at  the  center  of  the 
disc,  we  get  the  focussed  pressure  wave, 


p( 0)  ~  kvJo(0)  =  kv. 


M-k) 


&. 


Figure  1:  Plot  of  first  few  Bessel  Functions  Jn(z). 


Next  Steps:  Here  are  some  ideas  on  how  to  proceed  from  here  on: 
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•  Clearly  some  follow-thru  is  needed.  Short  computations  with  this  Bessel  solution  would  add  alot  to 
this  little  write-up. 

•  There  is  a  difficulty  in  the  cases  when  J\  (kR)  =  0  do  not  allow  a  solution  to  the  problem.  Probably 
insisting  that  a  >  0  will  remedy  this  but  the  analysis  needs  to  be  done.  Perhaps  a  is  small  enough, 
after  nondimensionalization,  to  develop  a  perturbation  solution  using  the  Bessel  solution  above. 


3.  Inverse  Problem: 

The  computation  of  various  constant  or  function  coefficients  in  ordinary  or  partial  differential  equation 
problems  is  commonly  referred  to  as  an  inverse  problem.  Since  the  terms  involving  the  unknowns  often 
do  not  ’’coerce”  the  coefficient  being  sought,  the  task  of  determining  it  is  ill-posed.  This  instability  may 
manifest  itself  as  non-uniqueness  or  extreme  sensitivity  to  noise/errors  in  the  given  data. 

There  are  many  approaches  to  solving  inverse  problems.  The  book  [G]  gives  an  overview  of  many  of  the 
issues.  The  paper  [T]  gives  one  example  of  an  ”ad  hoc”  technique  to  solve  an  inverse  problem  while  [EKN] 
and  [EZ]  provide  an  analysis  of  the  standard  linearization  approach.  The  paper  [Mu]  provides  an  entirely 
different  algorithm. 

We  hope  to  use  the  ideas  in  [EKN]  and  [EZ]  or  [Mu]  and  [T]  to  develop  a  theoretical  analysis  for  a 
technique  that  is  different  from  the  Bayesian  approach  in  [HHMK] . 

It  seems  reasonable  to  consider  a  one-dimensional  situation  (Proof  of  concept)  at  first.  We  should  aim  to 
develop  both  a  computer  solution  and  theory  about  the  problem.  In  one-dimension,  if  a  and  p  are  assumed 
piecewise  constant,  as  in  [HHMK],  then  a  solution  formula  for  p  =  p(x)  seems  possible  The  problem  defined 
by  (5)  and  (6)  is  primarily  steady  state.  However,  since  the  vn  may  depend  on  £,  apparently  they  ’’pulse” 
the  ultrasound,  the  solution  p,  in  fact,  could  change  its  steady  profile  from  pulse  to  pulse. 

We  think  it  is  of  value  to  investigate  the  qualitative  solution  behavior  using  either  numerics  or  closed 
form  solutions  probably  assuming  that  /?,  ft,  a  and  p  are  constant,  or  at  best,  piecewise  constant.  From  this 
we  will  have  a  clearer  notion  of  solution  behavior  and  dynamics. 

One-Dimensional  Elliptic  Problem  with  Piecewise  Constant  Parameter:  We  develop  and  analyze  a 
method  for  the  following  inverse  problem:  Find  function  [3  =  /3(x)  that  satisfies  the  boundary  value  problem 
(BVP) 

—U"  +  f3U  =  /  with  boundary  conditions  U( 0)  =  U(  1)  =  0. 

Suppose  there  is  a  piecewise  linear  data  function  Ud  defined  on  a  partition  of  Q  =  [0, 1]  with  subintervals 
of  uniform  width  hp  and  that  there  is  a  true  solution  of  the  BVP  U  =  U(x)  associated  with  this  [3 .  We  will 
assume  that  the  data  is  close  to  the  true  solution;  there  exists  0  <  e  <<  1  so  that 


\\U-UD II  <e. 


(7) 


We  further  assume  that 


f3(x)  = 


(3o  f°r  0  <  x  <  m, 

(3\  for  m  <  x  <  1, 


0  <  m  <  1. 


We  create  smooth  cutoff  functions  cjo  E  Co°(0,m)  and  E  Co°(ra,  1).  It  is  reasonable  to  require  that 
derivatives  of  the  uVs  are  0(1)  since  they  ’’live”  on  0(1)  intervals.  We  compute  an  approximation  f3A  using 
integration-by-parts  twice,  the  fact  that  cjo  and  are  smooth  and  have  support  inside  [0,  m]  and  [m,  1], 
respectively,  leads  to 


aA  _  In  M>  dx  +  fa  UDuo  dx  vA  _  In  M  dx  +  fa  u^'[  dx 

^  fQ  Uduiq  dx  Pl  fQ  UDu>\  dx 
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We  now  will  find  error  bounds  for  the  quantities  |/?o  —  (3q\  and  \(3i  —  /3f\  in  terms  of  hn  and  e.  It  will  be 
necessary  to  require  that  there  exist  positive  constants  po  and  pi  so  that 

|  /  Uljo  dx |  >  po  >  0  and  |  /  Uu\  dx\  >  pi  >  0. 

Jq  Jq 

We  are  pursuing  an  a  priori  style  result  and  so,  we  will,  implicitly  assume  that  norms  of  U  and  /  as  well  as 
bounds  of  /3q  and  f3\  are  all  0(1). 

Lemma  1:  There  exists  a  constant  C  independent  of  e  and  hn  such  that 


Ud^o  dx |  >  po  —  Ce 


and 


C/d  a; i  dx  |  >  pi  —  Ce. 


Proof:  Use  true  solution  U  as  an  intermediate  function; 


Po  <  |  /  Ux>0  dx\  <  I  /  (C/  -  Ud)uo  dx |  +  |  /  Udcjo  dx| 

JO  JO  JO 

Now,  using  the  estimate  7  we  can  finish  the  proof  for  the  Oth  term.  The  estimate  for  the  other  is  similar. 

/// 

Thus,  lemma  1  guarantees  that  /3q  and  [3 f  can  be  computed  since  the  denominators  in  their  definitions 
cannot  be  zero. 


Theorem:  There  exists  a  constant  C  independent  of  e  and  hp  such  that 

max{|/?0  ~  Po\,  \Pi  ~  Pi  \}  <Ce. 

Proof:  We  show  this  just  for  the  Oth  components.  The  other  term  is  similar.  We  start  by  subtracting  the 
key  integrated  equations  for  us  and  U  that  arose  from  their  respective  definitions; 

(  [  Uduj0  dx)  $-([  Ulo0  dx)  0o=  f  (UD  -  U)u'0'  dx 
JO  J o  /o 

or 

(  [  Udlo0  dx)(0$  -0o)=  f  (Ud  -  uy0'  dx  +  [  (UD  -  U)u0  dxpo ■ 

JO  JO  JO 

We  can  now  finish  our  proof  using  our  previous  lemma  and  assumption  that  we  are  pursuing  an  a  priori 
style  analysis.  /// 

Kristen’s  Multi-Parameter  Study:  The  above  has  been  extended  to  the  inverse  problem  below  where 
the  unknown  parameters  k  and  f3  are  sought  in  the  BVP 

-kU"  +pU  =  f  with  U{  0)  =  t/(  1)  =  0 

given  /  =  f(x)  and  data  for  U  =  Ud  where  /3  =  0(1000)  and  k  =  0(1).  Both  a  theoretical  analysis  and 
computations  with  finite  difference  schemes  has  been  done  to  show  the  veracity  of  the  approach. 

The  difference  in  sizes  for  n  and  /3  can  lead  to  subtractive  cancellation  issues  (due  to  data  approximation 
errors  (not  floating  point))  in  the  case  when  U  and  U"  are  0(1)  and  /  is  0(/3).  To  be  more  precise,  our 
direct  residual  technique  involves  creating  equations  for  the  parameters  by  multiplying  the  DE  by  functions 
l o  and  integrating  by  parts  over  U  =  [0, 1]; 

_  Snfu  d^^^Uixdx 
fQUu"  dx 
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So,  specifically,  the  subtractive  cancellation  arises  from  the  difference  calculation  of  the  two  large  but  close 
(in  value)  terms  in  the  numerator.  Choosing 


w(x)  =  e(l  -  e-x/e)(l  -  e-Q-*)h) 

(0  <  e  <<  1)  could  help  since 

[  Ulo  dx  =  O(e),  f  fix  dx  =  0(/3e)  and  [  U"ixdx  =  f  Uu"  dx  =  0(1) 

JQ,  Jo,  Jo,  JQ, 

so  the  calculation  of  / 3  involves  the  subtraction  of  smaller  numbers.  If  this  is  unsuccessful  the  other  possibility 
is  to  investigate  modifying  /  so  that  the  second  derivative  term  is  large  by  itself. 

Another  alternative  would  involve  approximating  the  difference  with  the  0(/3)  terms  leading  to  k  con¬ 
siderably  more  accurately. 

Full  Parameter  Set  Study  on  Reduced  Dimension  Problem:  The  next  step  would  be  to  add  the 

Helmholtz  equation  and  study  the  inverse  problems  with  unknown  parameters  o,  /?,  k  and  p.  The  most 
interesting,  and,  perhaps  realistic  version,  would  only  assume  there  was  data  for  the  temperature.  The  two 
equations  are  connected  thru  the  source  term  which  may  lead  to  further  ill-posedness  (see  page  64-65  in 
[G]).  This  analysis/development  would  involve  incrementally  more  complicated  problems,  probably  starting 
with  constant  coefficients  and  using  known  solutions  (Mainly  for  Helmholtz)  to  finite  difference  schemes  with 
meshes  that  fit  the  piecewise  constant  partitions  exactly  -  but  modest  frequencies  to,  finally,  more  realistic 
versions  using  the  enrichment  techniques  described  in  section  4. 

4.  Helmholtz  Equation  —  High  Frequencies  and  Rough  Coefficients: 

As  in  [HHMK]  we  use  the  Helmholtz  equation  to  model  the  pressure  instead  of  the  KZK  [HMB] ,  full  second- 
order  wave  [SZJ],  Burgers  or  Westervelt  equations.  (Investigation  of  these  other  models  in  this  MRI/HIFU 
context  might  provide  a  new  research  direction.). 

Challenges  &  Existing  Techniques:  There  are  several  major  challenges  associated  with  solving  the 
Helmholtz  equation  in  the  HIFU  context.  The  second  order  term  ” competes”  with  the  reaction  term  making 
standard  stability /energy  analyses  more  difficult.  The  absorption  term  (a  >  0)  alleviates  this  to  some  ex¬ 
tent.  At  the  same  time,  the  transducers  are  operated  in  at  high  frequencies;  for  accurate  resolution  in  finite 
difference  or  element  approaches,  very  fine  mesh  refinements  are  required  (See  [IH1]  and  [IH2]  for  more  on 
these  challenges). 

The  relatively  current  work  [MW]  uses  a  Treffitz  type  method.  A  discontinuous  finite  element-type  space 
is  developed  using  known  solutions  to  the  Helmholtz  equation  on  each  element  domain  -  triangles  usually  - 
in  a  least  squares  variational  form.  This  strategy  is  used  in  [HHMK];  the  meshes  are  fitted  to  the  different 
large  subdomains  (Skin,  tissue,  tumor  and  air).  Thus,  on  each  subdomain  a  constant  coefficient  Helmholtz 
equation  is  being  solved. 

Enrichment  with  XFEM  and/or  Pseudo-Derivative  Mesh  Free  Schemes:  Imbedding  appropriate 
singular  functions  into  approximation  spaces  is  a  common  idea  that  we  will  examine  to  handle  both  the 
material  discontinuity  and  high  frequency  challenges  (Although  the  approach  in  [MW]  handles  this  via  the 
introduction  of  Helmholtz  solutions  directly  into  the  local  element  functions).  We  are  interested  in  XFEM 
(see  [GWB],  [SVC]  or  [VSC])  as  well  as  meshfree  (see  [BM],  [FO],  [KLYBL]  or  [LY])  approaches. 

5.  Other  Directions: 
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Below  are  a  few  other  research  paths  we  have  pondered. 

Nonlinear  Inverse  Problem:  We  might  extend  the  current  work  by  looking  for  conductivity  and  perfusion 
as  functions  of  temperature;  that  is ,  k  =  k(T)  and  (3  =  /3(T).  (The  use  of  Grobner  bases  may  be  relevant, 
hard/interesting?  Need  to  look  into  work  by  Marisa  Eisenberg,  Hans  Othmer  and  Brandy  Stigler.). 

Optimization  of  Tranducer  Waves:  The  authors  of  [HHMK]  are  interested  in  the  transducer  array 
operation.  They  provide  an  optimization  algorithm  that  seeks  to  arrange  the  ultrasound  sources  to  provide 
focused  heating  but,  at  the  same  time,  avoid  overheating  -  to  minimize  tissue  damage. 

Sensitivity:  We  will  eventually  need  to  estimate  sensitivities  -  either  theoretically  or  computationally.  If  n 
is  an  unknown  then  we  would  be  trying  to  find  dn / 3U d  • 
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L 
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Ct 
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Ta 
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LOb 

Blood  Density 
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f 
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K, 
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0.6  W  K-1  m_1 

Un 

Transducer  Velocity 

1000  ms~x 
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Table  2  —  Estimates  of  Parameters  from  [HHMK 


(Page  1020): 


Material 

a  (dB/(MHz  cm)) 

c  (ms  -1) 

p  (kg  m~3) 

3? 

1 

1 

CT  (Jkg-'K-1) 

(jJb  (kg  m  3s  -1) 

Air 

0 

1500 

1000 

0.60 

4190 

0 

Skin 

12 

1610 

1200 

0.50 

3770 

1.0 

Tissue 

5 

1485 

1020 

0.50 

3550 

0.7 

Tumor 

5 

1547 

1050 

0.65 

3770 

2.3 

Appendix  I  -  Solutions  to  the  Linear  Pressure  Equation. 

If  we  consider  the  special  case  where  density  p  is  identically  equal  to  1  in  one-dimension  on  [0,  L\  we  can 
make  a  beginning  assessment  of  the  form  of  the  solution  p(x).  Here  we  wish  to  solve 

p"  +  k2p  =  0  p'(0)  =  iPo  (k  =  f  +  ia,  a  >  0.) 

If  the  boundary  condition  at  the  right  ”  end”  point  is  either  absorbing  as  in  [HHMK] , 

p'(L)  -  ikp(L)  =  0 


or  requires  p  — >  0  as  x  — »  oo,  then  we  have 

„  pnf>ifxp-ax 

p(x)  =  Aeifxe~ax  or  p(x)  =  —£ - ; - . 

J  +  lOL 

Note,  in  particular,  in  this  case  the  high  frequency  oscillation  is  not  present  in  the  heat  source  term  Q  =  Q(x) 
since  the  absolute  value  of  pressure  becomes: 


We  can  then  choose  Pq  =  Aofeax  so 


\p(x)\ 


P0e~aa 

\//2  +  a2 


\p{x)\  =  A0 


e-a(x-x  0) 
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If,  instead,  we  impose  p'(L)  =  0  at  the  right  end  point  we  find 


p(x) 


iPocos(k(x  —  L )) 
ksm(kL) 


and 


\p(x)\ 

MR-Guided  HIFU  System: 


Po  /  cos(2  fx)  +  cosh(2<xr) 
yJJ 2  +  a2  \j  cos(2  fL)  —  cosh(2 aL) 
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Gj 


Figure  2:  A  feedback  control  system  for  preclinical  thermal  therapy  was  established  using  MRguided  HIFU  in  the 
CCHMC  IRC.  The  feedback  control  loop  is  formed  using  a  3.0  Tesla  32-channel  Philips  Achieva  MR  scanner  (Philips 
HealthCare,  Best,  the  Netherlands)  and  a  Philips  small-animal  HIFU  system  (Philips  HealthCare,  Vantaa,  Finland). 
The  current  system  uses  a  single  channel  RF  coil  for  animal  imaging  and  an  eight-channel  3  MHz  sectored  transducer 
(IMASONIC,  Voray  sur  I’Ognon,  France)  for  acoustic  transmission.  The  development  of  a  multichannel  RF  coil  array 
for  high-speed  MR  imaging  of  mice  is  ongoing. 
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